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SUPPLEMENTAL INFORMATION 


Figure S1. Waveform and spectrogram of a Bornean orangutan long call, related to 


STAR methods. Rectangles indicate the selections bounding individual pulses. 


Selections were manually drawn on spectrograms by observers using Raven Pro 1.6 


software (K. Lisa Yang Center for Conservation Bioacoustics, 2019). Call duration was 


measured as the interval between the start of the first selection to the end of the last 


selection. 
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Figure S2. Distribution of focal observations by male and month. Stacked-bar plot 
depicting the number and distribution of full-day behavioral observations, disaggregated 


by male ID. 
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Table S1. Model outputs examining immediate (I) and persistent (P) effects of smoke on long calls, related to Figures 5 and 6. 
Mixed model outputs for biphonation (biph 0-1 = presence, % = relative duration), deterministic chaos (chaos 0-1 = presence, % = 
relative duration), fundamental frequency (f0), jitter, harmonics to noise ratio (HNR), shimmer (shim), daily call rate (N call/day), and 
call duration; all models included male ID as a random intercept (i.e., Variable ~ Period + (1 | Male.ID)). GLMMs used either binomial 
(bin) or poisson (pois) distributions. Estimates of the beta coefficients (Est.) for the effect of smoke period, their standard error (SE), t- 
statistic, and p-values for both 2-tailed (2-t) and 1-tailed (1-t) tests. Marginal (marg) R’ represents the variance explained by fixed 
factors and conditional (cond) represents the variance explained by both fixed and random factors. For call duration, we report model 
outputs from two models: one including only period and one including both period and call context (spont) as fixed factors. Because 
our 3-male model for persistent f0 (3-male) produced a ‘singular fit’ warning (indicated overfit), we removed one male (Tomi) from 


the analysis and ran the model again (2-male). All other models converged as expected. 


Period 


I/P_ Variable Data Model Type N Est. SE zit P(2-t) P(i-t) | marg. cond. 
I Biph. (0-1) nlp_PreSmoke GLMM (bin) 136 1.010 0.392 2.576 0.010 0.005 | 0.049 0.178 
I Biph. (%) nlp_PreSmoke LMM 136 | -0.004 0.051 -0.079 0.937 1.875 | 0.000 0.126 
I Chaos (0-1) nlp_PreSmoke GLMM (bin) 136 0.185 0.365 0.506 0.613 0.306 | 0.002 0.002 
I Chaos (%) nlp_PreSmoke LMM 135 0.031 0.046 0.677 0.500 0.250 | 0.003 0.003 
I f0 (call avg) voice PreSmoke LMM 108 | -13.084 6.228 -2.101 0.038 0.019 | 0.040 0.042 
I jitter (call avg) voice PreSmoke = LMM (log) 108 0.006 0.003 2.168 0.032 0.016 | 0.026 0.365 
I HNR (call avg) = voice PreSmoke LMM 108 | -1.115 0.663 -1.683 0.095 0.048 | 0.019 0.315 
I shim (call avg) voice PreSmoke LMM 108 0.015 0.010 1.478 0.143 0.071 | 0.011 0.341 
I N call/day N.LC_PreSmoke |= GLMM (pois) 56 -0.508 0.143 -3.540 0.000 0.000 | 0.104 0.623 
I Call duration durat_PreSmoke LMM 119 | -3.033 3.266 -0.929 0.355 0.178 | 0.007 0.109 
I Call duration durat_PreSmoke | LMM (spont) 119 | -2.762 3.151 -0.876 _ 0.383 0.191 | 0.078 0.166 
P ___Biph. (0-1) nlp_PrePost GLMM (bin) 79 -1.845 0.631 -2.922 0.003 0.002 | 0.112 0.143 
P __Biph. (%) nlp_PrePost LMM 79 0.124 0.0706 1.75 0.084 0.168 | 0.037 0.084 
P Chaos (0-1) nlp_PrePost GLMM (bin) 19 0.642 0.625 1.028 0.304 0.608 | 0.016 0.016 
P Chaos (%) nlp_PrePost LMM 79 0.052 0.057 0.912 0.365 0.729 | 0.011 0.011 
P — f0 (call avg) voice_PrePost LMM (3 male) 58 10.472 10.594 0.988 0.327 0.164 | 0.017 0.017 
P f0 (call avg) voice_PrePost LMM (2 male) 49 12.735 10.871 1.171 0.247 0.124 | 0.028 0.030 
P jitter (call avg) voice_PrePost LMM (log) 58 -0.012 0.005 -2.473 0.017 0.008 | 0.069 0.227 
P  HNR(callavg) voice _PrePost LMM 58 3.593 1.199 2.997 0.004 0.002 | 0.117 0.261 
P shim (call avg) voice_PrePost LMM 58 -0.042 0.017 -2.532 0.014 0.007 | 0.047 0.092 
PN call/day N.LC_PrePost GLMM (pois) 45 2.311 0.324 7.145 0.000 0.000 | 0.745 0.866 
P Call duration durat_PrePost LMM 64 11.254 5.924 1.900 0.062 0.031 | 0.047 0.189 
P Call duration durat_PrePost LMM (spont) 64 11.060 5.784 1.912 0.061 0.030 | 0.099 0.235 


Table S2. Sample sizes for NLP analysis, related to STAR methods and Figure 3. 


Number of calls and pulses manually annotated for biphonation and deterministic chaos, 


disaggregated by male ID and study period. 


Pre-Smoke Smoke Post-Smoke 
Male ID Calls Pulses Calls Pulses Calls Pulses 
Niko 8 20 6 19 1 5 
Otto 8 19 6 12 - - 
Tomi 8 24 7 13 1 6 
Wodan 9 15 11 14 4 9 


Table S3. Descriptive statistics for environmental and contextual variables, related 
to STAR methods. Mean and range of maximum temperature (°C), day rainfall 


accumulation (mm), and spontaneous calls (proportion of all calls) disaggregated by 


study period. 
Pre-Smoke Smoke Post-Smoke 
Variable Mean Range | Mean Range | Mean Range 
Max temp 31.21 26-34 | 28.75 27-31 | 31.21 26-34 
Day rain 118 0-82 0 - 1.41 0-7.8 
Spontaneous 0.75 - 0.76 - 0.78 - 


Table S4. Model outputs examining environmental and contextual effects on long 


call features, related to STAR methods. Estimates of the beta coefficients (Est.) for the 


effect of each variable, their standard error (SE), z- or t-statistic (z/t), and p-value (p) on 


daily number of long calls (N call) and mean duration, FO, jitter, hnr, and shimmer. Each 


model included one environmental or contextual predictor variable and smoke period as 


fixed factors plus male ID as a random intercept (i.e., Response. Variable ~ Period + 


Predictor. Variable + (1 | Male.ID)). Models were constructed separately for immediate 


(i.e., during smoke period = I) and persistent (1.e., during post-smoke period = P) changes 


relative to the pre-smoke period. All models examining FO (except immediate-rain) 


produced a ‘singular fit’ warning (indicating overfit). 


N call duration FO 
Factor | /P | Est. SE zit p Est. SE zi/t p Est. SE  z/t p 
Temp I | -0.02 0.08 -0.25 0.800) -0.34 1.73 -0.20 0.846 | 7.83 3.87 2.02 0.046 
Rain I | -0.03 0.03 -0.85 0.395] -0.94 1.40 -0.68 0.501] 0.91 2.59 0.35 0.725 
Spont I - - - - -11.58 3.67 -3.15 0.002 | 11.76 7.44 1.58 0.117 
Temp P |} -0.01 0.07 -0.11 0.911 | -068 2.00 -0.34 0.735 | 2.73 3.44 0.79 0.431 
Rain P |} -0.01 0.03 -0.47 0.638] -0.85 1.05 -0.81 0.423 | -0.53 1.92 -0.28 0.781 
Spont | P - - - - -9.52 4.74 -2.01 0.049 | 15.68 9.21 1.70 0.094 

jitter hnr shimmer 
Factor | V/P| Est. SE  z/t p Est. SE zi/t p Est. SE  z/t p 
Temp I | 0.00 0.00 0.81 0.421) -0.63 0.49 -1.29 0.199} 0.00 0.00 0.90 0.370 
Rain I | 0.00 0.00 -0.16 0.877) -0.02 0.28 -0.06 0.954} 0.00 0.00 -0.72 0.472 
Spont I | 0.00 0.00 0.02 0.982 | -0.28 0.80 -0.35 0.725 | 0.01 0.01 0.79 0.430 
Temp P | 0.00 0.00 0.98 0.331] -046 0.39 -1.18 0.243) 0.01 0.00 1.73 0.089 
Rain P | 0.00 0.00 -1.01 0.315] 0.11 0.22 048 0.631] 0.00 0.00 -0.94 0.350 
Spont | P | 0.00 0.00 0.08 0.936} -0.30 1.10 -0.27 0.789) 0.01 0.01 0.59 0.559 


